Exploration of Damage Identification Method for a Large-Span Timber Lattice Shell Structure in Taiyuan Botanical Garden based on Structural Health Monitoring

Large-span spatial lattice structures generally have characteristics such as incomplete modal information, high modal density, and high degrees of freedom. To address the problem of misjudgment in the damage detection of large-span spatial structures caused by these characteristics, this paper proposed a damage identification method based on time series models. Firstly, the order of the autoregressive moving average (ARMA) model was selected based on the Akaike information criterion (AIC). Then, the long autoregressive method was used to estimate the parameters of the ARMA model and extract the residual sequence of the autocorrelation part of the model. Furthermore, principal component analysis (PCA) was introduced to reduce the dimensionality of the model while retaining the characteristic values. Finally, the Mahalanobis distance (MD) was used to construct the damage sensitive feature (DSF). The dome of Taiyuan Botanical Garden in China is one of the largest non-triangular timber lattice shells worldwide. Relying on the structural health monitoring (SHM) project of this structure, this paper verified the effectiveness of the damage identification model through numerical simulation and determined the damage degree of the dome structure through SHM measurement data. The results demonstrated that the proposed damage identification method can effectively identify the damage of large-span timber lattice structures, locate the damage position, and estimate the degree of damage. The constructed DSF had relatively strong robustness to small damage and environmental noise and has practical application value for SHM in engineering.


Introduction
During long-term service, structures may suffer damage due to material aging, corrosion, prolonged loading, and natural disasters such as fire and earthquakes [1]. For major engineering structures, such as large-span bridges, high-rise buildings, large-span spatial lattice structures, large-scale water conservancy structures, and large-scale offshore platforms, failure can cause significant economic damage and human casualties. Therefore, SHM is of great significance in civil engineering and related fields [2].
Damage identification and localization are crucial processes in SHM since damage alters the stiffness, mass, or damping of structures. Consequently, stiffness reduction caused by cracking or loosening of connections will lead to changes in structural vibration modes, resulting in alterations in its dynamic response [3,4]. However, obtaining complete and accurate high-order modal information through theoretical analysis for actual structures, particularly for large-span spatial structures with numerous degrees of freedom, is challenging due to the complex vibration data involved. In addition, in practical engineering, it is almost impossible to obtain input and output signals by artificially exciting the structure and accurately identifying the structural modal information.
structures, which does not require modal parameter identification and excitation information, avoids the dependence of damage identification results on the accuracy of structural finite element (FE) models, and has high damage identification accuracy.
This article introduced PCA on the basis of existing research to extract features from the residual sequence. The advantage of incorporating PCA lies in its ability to reduce the dimensionality of the data while retaining the most informative features [17]. This not only simplifies the subsequent analysis but also helps to identify the most significant damage indicators. In addition, the MD was used in constructing DSF because it takes into account the differences and correlations in the variability of each observed variable. In addition, it can effectively eliminate the influence of scale on different measurement values.
Furthermore, existing studies were limited to theoretical research and simple experimental verification, lacking engineering applications and practices for large and complex structures. Based on the SHM project for a timber lattice shell dome in Taiyuan Botanical Garden, this paper proposes a damage identification method based on time series analysis. The effectiveness of the damage identification model was verified using numerical simulation, and the degree of damage to the lattice shell structure was identified using measured data from SHM. This paper combined theoretical research with engineering applications to verify the effectiveness and feasibility of applying the time-series-analysisbased damage identification method for practical engineering.

Research Background
This study took the timber lattice shell of Taiyuan Botanical Garden as the research background, which is located in Taiyuan, Shanxi, China. It is one of the largest non-triangular timber lattice structures worldwide. It was designed by Delugan Meissl Associated Architects (DMMA) from Austria and is now a landmark building in Taiyuan. The dome, from above, is shaped like a shell and glazed with double-curved panes of glass. The span of the dome is 89.5 m, with a span-to-height ratio of 3.0 and a projected area of approximately 6000 m 2 . The dome consists of double-curved laminated timber beams, which are arranged in two or three intersecting layers. The beams in the intersection joint area are connected by steel pins and self-tapping screws. The splicing connection of the beam adopts the half-lap joint form and is connected by self-tapping screws. To improve the overall rigidity of the structure, a bidirectional prestressed steel cable net was installed inside it, arranged diagonally to the timber beams. The steel support of the dome is connected to the timber beam using anti-shear steel plates and self-tapping screws and connected to the concrete foundation using chemical anchor bolts. The appearance and structural details of the dome are shown in Figure 1. For more details of the dome, please refer to reference [18].

Methodology
This section provides an overview of the fundamental principles of utilizing the ARMA model for time series analysis. It also elucidates the techniques for determining the order and estimating the parameters of the ARMA model using AIC and LAR. Furthermore, this section offers a theoretical exposition of the principal component analysis method introduced and defines a DSF construction method based on MD.

Time Series Analysis Modeling
Time series analysis is a statistical method for dynamic data processing, which is based on the theory of stochastic processes and mathematical statistics. It uses parameter models to process ordered random sampling data and studies the statistical laws followed by data sequences, in order to perform system identification [19,20]. Any structure can be regarded as a mechanical system composed of stiffness, mass, and damping matrices, which can be described by a motion differential equation and can be transformed into a difference equation in the discrete time domain. This difference equation has the same expression form as the time series model. Therefore, the equivalent relationship between the time series model and the structural motion differential equation can be used to determine the state of the structure and perform damage identification of the system.

Basic Theory of ARMA Model
The ARMA model is a statistical model commonly used in time series analysis. It takes past observations of the time series as independent variables and current observations as the dependent variables to describe the randomness and autocorrelation of the time series data. The ARMA model consists of two parts: the AR model and the moving average (MA) model, and has the characteristics of both models. The AR or MA model is a special case of the ARMA model, so the ARMA model is more general for system response and is the most commonly used model in time series analysis.
For a multi-degree-of-freedom system, its vibration differential equation is: where M, C, and K represent the mass matrix, damping matrix, and stiffness matrix of the structure, respectively. ( ), ( ), and ( ) represent the acceleration vector, velocity vector, and displacement vector of the structure, respectively. ( ) is the excitation force applied to the system. Equation (1) is equivalent to a 2n-order non-homogeneous differential equation sys-

Methodology
This section provides an overview of the fundamental principles of utilizing the ARMA model for time series analysis. It also elucidates the techniques for determining the order and estimating the parameters of the ARMA model using AIC and LAR. Furthermore, this section offers a theoretical exposition of the principal component analysis method introduced and defines a DSF construction method based on MD.

Time Series Analysis Modeling
Time series analysis is a statistical method for dynamic data processing, which is based on the theory of stochastic processes and mathematical statistics. It uses parameter models to process ordered random sampling data and studies the statistical laws followed by data sequences, in order to perform system identification [19,20]. Any structure can be regarded as a mechanical system composed of stiffness, mass, and damping matrices, which can be described by a motion differential equation and can be transformed into a difference equation in the discrete time domain. This difference equation has the same expression form as the time series model. Therefore, the equivalent relationship between the time series model and the structural motion differential equation can be used to determine the state of the structure and perform damage identification of the system.

Basic Theory of ARMA Model
The ARMA model is a statistical model commonly used in time series analysis. It takes past observations of the time series as independent variables and current observations as the dependent variables to describe the randomness and autocorrelation of the time series data. The ARMA model consists of two parts: the AR model and the moving average (MA) model, and has the characteristics of both models. The AR or MA model is a special case of the ARMA model, so the ARMA model is more general for system response and is the most commonly used model in time series analysis.
For a multi-degree-of-freedom system, its vibration differential equation is: where M, C, and K represent the mass matrix, damping matrix, and stiffness matrix of the structure, respectively. ..
x(t), and x(t) represent the acceleration vector, velocity vector, and displacement vector of the structure, respectively. F(t) is the excitation force applied to the system. Equation (1) is equivalent to a 2n-order non-homogeneous differential equation system. Under the action of a single excitation f (t), the vibration differential equation is expressed as follows: ϕ 2n x (2n) + ϕ 2n−1 x (2n−1) + · · · + ϕ 1 . x Discretizing Equation (2) and using a white noise sequence a t ∼ ND 0, σ 2 a ) as the system input yields a single degree of freedom ARMA (p,q) model [21]: In the equation, ϕ i represents the i-th order AR coefficient, θ j represents the j-th order MA coefficient, and p and q are the AR and MA orders, respectively. By introducing a backward shift operator B, defined as B k x t = x t−k [5], Equation (3) can be transformed into: According to Equation (4), the ARMA model describes a system with a transfer function of θ(B)/ϕ(B), where ϕ(B) represents the inherent characteristics of the system and θ(B) represents the relationship between the system and the external environment. As the input signals of large civil engineering structures are often difficult to test or measure, the performance and behavior of the structure can only be inferred through the analysis of output signals. The ARMA model can be established by analyzing the output signal of the structure, without considering the specific information of the input signals, and only using the white noise sequence {a t } as input to establish an analysis model. This makes the ARMA model increasingly widely used in the field of civil engineering. In the analysis of time series models, the AR, MA, and ARMA models have their own characteristics, which can be used to preliminarily determine the appropriate time series model. The characteristics of the models are shown in Table 1. When applying the ARMA model for time series analysis, it is necessary to determine the order of the ARMA model, and the quality of the order setting has a significant impact on the parameter identification results. If the order is too high, the model will be too complex, resulting in overfitting and unnecessary calculations; if the order is too low, the model will be too simple to capture important features of the time series, leading to poor model fitting. Currently, one of the mostly used methods for determining the order of the ARMA model is AIC [22]. The function of AIC is defined as follows: In Equation (5), N represents the length of the time series x t , σ 2 a represents the variance of the ARMA model residuals, and p and q represent the AR and MA order, respectively. As the sum of model orders (p + q) increases, ln σ 2 a will decrease and 2(p + q)/N will increase. Therefore, under the given parameter estimation method, the order with the minimum AIC(p + q) value should be chosen.

Estimation of ARMA Parameter
This article used the long autoregressive (LAR) method to estimate the parameters of the ARMA model. Compared with other methods for parameter estimation, the LAR method can transform nonlinear regression problems into linear regression problems, which can simplify the calculation process and improve efficiency. According to the theory of LAR method, the ARMA model and the AR model are equivalent mathematical models. Therefore, the parameters of the AR model can be estimated first, and then the parameters of the ARMA model can be estimated based on the relationship of the transfer functions. The equivalent relationship between the transfer functions of AR and ARMA models is as follows: In Equation (6), the left side represents the transfer function of an AR(p) model, and the inverse function I i is equal to the parameter ϕ i of the AR model. The right side represents the transfer function of an ARMA(n,m) model. By moving the same power coefficients of the shift operator B in Equation (6), it can be deduced that According to Equation (7), the solution for ϕ i (i = 1, 2, . . . , n) can be obtained as follows: For the last term in Equation (7), let k = n + 1, n + 2, . . . , n + m, where n + m = p. Then, I i (i = n + 1, n + 2, . . . , n + m) can be expressed as Both Equations (8) and (9) are linear equation sets about θ j , so θ j can be solved first according to Equation (9), and then ϕ i can be solved according to Equation (8).
Based on the theoretical research of Ljung [23], this paper adopted the order scheme of p = 2m, q = 2m − 1 to establish the ARMA model of a multi-degree-of-freedom system, which is also a widely used model scheme in engineering [24].

Principal Component Analysis
There is a large amount of environmental noise in the on-site measurement data of a structure, which greatly reduces the efficiency and reliability of damage monitoring. Therefore, it is necessary to conduct feature analysis on the data to remove the noise influence and amplify the impact of structural damage on data changes. In this paper, PCA was adopted for feature extraction, which can filter out noise and redundancy in the original high-dimensional feature space data and transform it into interpretable data in a low-dimensional feature space. It can eliminate confusing data while preserving the main information [25].
For a dataset X ∈ R n×p with p variables and n samples, the i-th principal component of X can be represented as where u pi represents the characteristic vector. After standardizing the original data matrix X, the resulting matrix X* has a covariance matrix S that is equal to the correlation coefficient matrix: The eigenvalues of the covariance matrix of X* are denoted by λ 1 , λ 2 , . . . , λ p , and The number of principal components n is determined by the variance accumulative contribution rate (ACR) If the ACR of the first n principal components is not less than 85%, it indicates that the first n principal components have incorporated most of the original data information [26].

Damage Sensitive Feature
DSF refers to the characteristic or indicators used to identify and locate potential structural damage, which can reflect the internal state changes of the structure or system. It requires the establishment of an ARMA model to analyze the signal, and the parameters of the model include the inherent characteristics of the system. By extracting the carrier containing structural features, it can be determined whether the structure is damaged. Because the AR coefficients reflect the dynamic characteristics of the system, they can be used to construct the DSF.

Mahalanobis Distance
The principle of the structural damage detection method based on statistical pattern recognition is to compare two sets of model parameters, one set is the model feature parameters of the structure in the healthy state, and the other set is the feature parameters in the unknown (to be identified) state. In the damage identification method based on the ARMA model, the MD is usually used to determine the distance between the sample data and the reference data after feature extraction, and then evaluate the health status of the structure. MD calculates the distance between two multivariate sets by considering the correlation between them, which is an effective method for computing the similarity between two unknown-state sample sets. The expression of the MD between the sample set z and the population G is: where v G represents the mean vector of G, and Σ represents the covariance matrix of G.

Construction of DSF
ARMA models were established for the unknown-state sample set, reference sample set, and training sample set. The DSF was constructed by measuring the differences between the data parameters of the ARMA model under the unknown state and the healthy state. The formula is as follows: where u represents the unknown-state sample set, r represents the reference sample set of the healthy state, and T represents the training sample set of the healthy state. A DSF value close to 1 indicates that the structure is in a healthy state, while a DSF value greater than 1 indicates damage. Moreover, the severity of the damage is positively correlated with the magnitude of the DSF.

Results and Discussions
This section provides an explanation and discussion of the damage identification results using the FE models and SHM measured data. Firstly, a concise introduction is provided for the FE model, which includes the setup of five simulated DCs. The establishment of a time series model and the detailed results of damage identification are then elaborated. Finally, the obtained results are thoroughly explained and discussed. Furthermore, this section is based on two years of SHM data collected from the timber lattice shell dome of Taiyuan Botanical Garden. The proposed damage identification method from this article is applied to identify the damage in the dome, and the identification results are subsequently analyzed and discussed.

Damage Identification Based on Numerical Models
Structural damage in practical engineering exhibits characteristics such as randomness and finiteness, making it difficult to determine the effectiveness and availability of time series models through measured data. However, FE methods can be used to simulate the health and damage status of the structure, in order to establish and validate the ARMA model.
To verify the feasibility and applicability of the ARMA model for large-span timber lattice structures, this paper used Midas Gen, a general structural analysis software, to perform FE modeling and analysis for the timber lattice dome of Taiyuan Botanical Garden. The stiffness of individual or multiple components was reduced to simulate the structural damage conditions (DC). Gaussian white noise was used to excite the structure, and the response results of multiple measuring points (MP) were extracted to evaluate the damage identification effect of the ARMA model.

Finite Element Model
This paper used 3D line elements to analyze the spatial structure of the dome. The timber beams were modeled with beam elements, while the cable net was modeled with cable (tensile only) elements. In this study, the maximum main stress and maximum vertical deformation were employed as the convergence indicators. The analysis findings demonstrated that as the number of elements increased from 7732 to 334,817, both the maximum stress and maximum deformation exhibited a gradual convergence. The mesh convergence plot is shown in Figure 2. Notably, the maximum stress was more sensitive to the element number, thus serving as the primary controlling factor. Taking into account both computational accuracy and efficiency, the 134,372-element model was ultimately selected, which achieved a maximum stress convergence within a 5% range and satisfied the required computational accuracy.
The FE model took into account the fact that the three-layer timber beams are not coplanar, as shown in Figure 3. The timber beams are made of European spruce glued laminated timber GL28h, with a main beam section of 200 mm × 400 mm (width × height) and a secondary beam section of 200 mm × 300 mm (width × height). The components of the lattice structure are mainly subjected to axial forces, and glued laminated timber has relatively ideal material properties. According to material tests [27] and related studies [28][29][30], the FE material model adopted an orthotropic bilinear elastic-plastic model. The material test data and FE material model were shown in Figure 4, and the material properties are detailed in Table 2. The cable has a diameter of 26 mm, and the material used is austenitic steel 06Cr17Ni12Mo2 (316 stainless steel), with an elastic modulus of 160,000 MPa and a design pre-tension value of 40 kN. The FE model took into account the fact that the three-laye coplanar, as shown in Figure 3. The timber beams are made of laminated timber GL28h, with a main beam section of 200 mm × 4 and a secondary beam section of 200 mm × 300 mm (width × heig the lattice structure are mainly subjected to axial forces, and glue relatively ideal material properties. According to material tests [28][29][30], the FE material model adopted an orthotropic bilinear e material test data and FE material model were shown in Figure 4 erties are detailed in Table 2. The cable has a diameter of 26 mm, austenitic steel 06Cr17Ni12Mo2 (316 stainless steel), with an ela MPa and a design pre-tension value of 40 kN. The FE model took into account the fact that the three-layer timber beams are not coplanar, as shown in Figure 3. The timber beams are made of European spruce glued laminated timber GL28h, with a main beam section of 200 mm × 400 mm (width × height) and a secondary beam section of 200 mm × 300 mm (width × height). The components of the lattice structure are mainly subjected to axial forces, and glued laminated timber has relatively ideal material properties. According to material tests [27] and related studies [28][29][30], the FE material model adopted an orthotropic bilinear elastic-plastic model. The material test data and FE material model were shown in Figure 4, and the material properties are detailed in Table 2. The cable has a diameter of 26 mm, and the material used is austenitic steel 06Cr17Ni12Mo2 (316 stainless steel), with an elastic modulus of 160,000 MPa and a design pre-tension value of 40 kN.    Based on the measured data of SHM, simple modifications were attempted on th model. Initially, the acceleration signal during the early stage of structural service was cessed and analyzed. The time-domain data were then transformed into frequency-dom data using Fourier transform. To eliminate high-frequency noise, a low-pass filter of 5 was applied. The power spectral density (PSD) of the acceleration data was calculated, the dominant frequency with significant amplitude in the signal was extracted through p detection. This allowed for the determination of the natural frequency of the structure.
Subsequently, modifications were made to the FE model parameters. The el modulus and density of the material were adjusted, along with the magnitude of the ditional constant load (converted into mass) on the roof, and the connection stiffne the joints. These modifications aimed to align the natural vibration characteristics o FE model with the analysis results obtained from the measured data.

Damage Conditions
When establishing a time series model, the first step is to determine the sampling interval Δ and the sample length L for the continuous signal. When sampling a continu signal, a reasonable sampling time interval is Δ < 1/(2 ), and a reasonable sam length is > 1/Δ , where is the highest frequency in the interested vibra modes and Δ is the minimum difference between adjacent frequencies [32]. Given the complex natural vibration characteristics of the lattice shell structure, a order natural vibration modal analysis was conducted for the FE model. The results o modal analysis are presented in Figure 5 and Table 3. The analysis revealed that the na frequencies of the lower order modes are closely spaced, with a difference of only 8.92 between the first and 70th modes. By the 94th mode, 90% of the cumulative participa mass of the vibration modes was reached, indicating that the first 94 modes could cap the majority of the natural vibration characteristics of the structure. Therefore, the appro   Based on the measured data of SHM, simple modifications were attempted on the FE model. Initially, the acceleration signal during the early stage of structural service was processed and analyzed. The time-domain data were then transformed into frequencydomain data using Fourier transform. To eliminate high-frequency noise, a low-pass filter of 50 Hz was applied. The power spectral density (PSD) of the acceleration data was calculated, and the dominant frequency with significant amplitude in the signal was extracted through peak detection. This allowed for the determination of the natural frequency of the structure.
Subsequently, modifications were made to the FE model parameters. The elastic modulus and density of the material were adjusted, along with the magnitude of the additional constant load (converted into mass) on the roof, and the connection stiffness of the joints. These modifications aimed to align the natural vibration characteristics of the FE model with the analysis results obtained from the measured data.

Damage Conditions
When establishing a time series model, the first step is to determine the sampling time interval ∆t and the sample length L for the continuous signal. When sampling a continuous signal, a reasonable sampling time interval is ∆t < 1/(2 f max ), and a reasonable sample length is L > 1/∆ f min , where f max is the highest frequency in the interested vibration modes and ∆ f min is the minimum difference between adjacent frequencies [32].
Given the complex natural vibration characteristics of the lattice shell structure, a 100-order natural vibration modal analysis was conducted for the FE model. The results of the modal analysis are presented in Figure 5 and Table 3. The analysis revealed that the natural frequencies of the lower order modes are closely spaced, with a difference of only 8.92 Hz between the first and 70th modes. By the 94th mode, 90% of the cumulative participation mass of the vibration modes was reached, indicating that the first 94 modes could capture the majority of the natural vibration characteristics of the structure. Therefore, the appropriate sampling frequency and sample length were calculated based on the first 94 modes. The maximum frequency in the first 94 modes was approximately 49.85 Hz, and the minimum difference between adjacent frequencies was 0.0109 Hz. According to the sampling theories [32], the sampling time interval of this structure should be less than 0.0101 s, and the sample length should be greater than 91.75 s. Therefore, a sampling frequency of 100 Hz and a sample length of 120 s were adopted to establish a time series model. According to the basic principle of time series analysis, the signal input needs to meet the characteristics of stationary, normal, and zero-mean, and Gaussian white noise can just meet the above characteristics. Thus, Gaussian white noise generated by MATLAB was used as the environmental excitation to conduct the time history analysis on the FE model. sample length of 120 s were adopted to establish a time series model. According to the basic principle of time series analysis, the signal input needs to meet the characteristics of stationary, normal, and zero-mean, and Gaussian white noise can just meet the above characteristics. Thus, Gaussian white noise generated by MATLAB was used as the environmental excitation to conduct the time history analysis on the FE model.    In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4 In practical engineering, vibration sensors were used for SHM of the dome s Specifically, nine vibration MPs were set up at the joint positions, taking into acc shape and characteristics of the structure. From above, the MPs were distribute central point of the dome and in the directions of east, west, south, and north. In tical direction, the MPs were evenly distributed at the top, middle, and bottom structure. During the time history analysis on the FE models, acceleration respo were extracted from the nine MPs. The layout of the MPs was shown in Figure  FE analysis process, stiffness reduction was applied to local timber beam and st elements to simulate damage or stiffness degradation of the components. This p veloped five sets of DCs, including damage to single and multiple components. simulated DCs are detailed in Table 4.   Table 4. In this paper, the process of establishing a time series mo  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure Specifically, nine vibration MPs were set up at the joint positions, taking into account th shape and characteristics of the structure. From above, the MPs were distributed at th central point of the dome and in the directions of east, west, south, and north. In the ver tical direction, the MPs were evenly distributed at the top, middle, and bottom of th structure. During the time history analysis on the FE models, acceleration response dat were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In th FE analysis process, stiffness reduction was applied to local timber beam and steel cabl elements to simulate damage or stiffness degradation of the components. This paper de veloped five sets of DCs, including damage to single and multiple components. The fiv simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of t Specifically, nine vibration MPs were set up at the joint positions, takin shape and characteristics of the structure. From above, the MPs were central point of the dome and in the directions of east, west, south, and tical direction, the MPs were evenly distributed at the top, middle, structure. During the time history analysis on the FE models, accelera were extracted from the nine MPs. The layout of the MPs was shown FE analysis process, stiffness reduction was applied to local timber be elements to simulate damage or stiffness degradation of the compone veloped five sets of DCs, including damage to single and multiple com simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were Specifically, nine vibration MPs were set up at the jo shape and characteristics of the structure. From abo central point of the dome and in the directions of eas tical direction, the MPs were evenly distributed at structure. During the time history analysis on the FE were extracted from the nine MPs. The layout of the FE analysis process, stiffness reduction was applied elements to simulate damage or stiffness degradatio veloped five sets of DCs, including damage to single simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure Specifically, nine vibration MPs were set up at the joint positions, taking into account th shape and characteristics of the structure. From above, the MPs were distributed at th central point of the dome and in the directions of east, west, south, and north. In the ver tical direction, the MPs were evenly distributed at the top, middle, and bottom of th structure. During the time history analysis on the FE models, acceleration response dat were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In th FE analysis process, stiffness reduction was applied to local timber beam and steel cabl elements to simulate damage or stiffness degradation of the components. This paper de veloped five sets of DCs, including damage to single and multiple components. The fiv simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of t Specifically, nine vibration MPs were set up at the joint positions, takin shape and characteristics of the structure. From above, the MPs were central point of the dome and in the directions of east, west, south, and tical direction, the MPs were evenly distributed at the top, middle, structure. During the time history analysis on the FE models, accelera were extracted from the nine MPs. The layout of the MPs was shown FE analysis process, stiffness reduction was applied to local timber be elements to simulate damage or stiffness degradation of the compone veloped five sets of DCs, including damage to single and multiple com simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were Specifically, nine vibration MPs were set up at the jo shape and characteristics of the structure. From abo central point of the dome and in the directions of eas tical direction, the MPs were evenly distributed at structure. During the time history analysis on the FE were extracted from the nine MPs. The layout of the FE analysis process, stiffness reduction was applied elements to simulate damage or stiffness degradatio veloped five sets of DCs, including damage to single simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure Specifically, nine vibration MPs were set up at the joint positions, taking into account th shape and characteristics of the structure. From above, the MPs were distributed at th central point of the dome and in the directions of east, west, south, and north. In the ver tical direction, the MPs were evenly distributed at the top, middle, and bottom of th structure. During the time history analysis on the FE models, acceleration response dat were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In th FE analysis process, stiffness reduction was applied to local timber beam and steel cabl elements to simulate damage or stiffness degradation of the components. This paper de veloped five sets of DCs, including damage to single and multiple components. The fiv simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of t Specifically, nine vibration MPs were set up at the joint positions, takin shape and characteristics of the structure. From above, the MPs were central point of the dome and in the directions of east, west, south, and tical direction, the MPs were evenly distributed at the top, middle, structure. During the time history analysis on the FE models, accelera were extracted from the nine MPs. The layout of the MPs was shown FE analysis process, stiffness reduction was applied to local timber be elements to simulate damage or stiffness degradation of the compone veloped five sets of DCs, including damage to single and multiple com simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were Specifically, nine vibration MPs were set up at the jo shape and characteristics of the structure. From abo central point of the dome and in the directions of eas tical direction, the MPs were evenly distributed at structure. During the time history analysis on the FE were extracted from the nine MPs. The layout of the FE analysis process, stiffness reduction was applied elements to simulate damage or stiffness degradatio veloped five sets of DCs, including damage to single simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure Specifically, nine vibration MPs were set up at the joint positions, taking into account th shape and characteristics of the structure. From above, the MPs were distributed at th central point of the dome and in the directions of east, west, south, and north. In the ver tical direction, the MPs were evenly distributed at the top, middle, and bottom of th structure. During the time history analysis on the FE models, acceleration response dat were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In th FE analysis process, stiffness reduction was applied to local timber beam and steel cabl elements to simulate damage or stiffness degradation of the components. This paper de veloped five sets of DCs, including damage to single and multiple components. The fiv simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of t Specifically, nine vibration MPs were set up at the joint positions, takin shape and characteristics of the structure. From above, the MPs were central point of the dome and in the directions of east, west, south, and tical direction, the MPs were evenly distributed at the top, middle, structure. During the time history analysis on the FE models, accelera were extracted from the nine MPs. The layout of the MPs was shown FE analysis process, stiffness reduction was applied to local timber be elements to simulate damage or stiffness degradation of the compone veloped five sets of DCs, including damage to single and multiple com simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were Specifically, nine vibration MPs were set up at the jo shape and characteristics of the structure. From abo central point of the dome and in the directions of eas tical direction, the MPs were evenly distributed at structure. During the time history analysis on the FE were extracted from the nine MPs. The layout of the FE analysis process, stiffness reduction was applied elements to simulate damage or stiffness degradatio veloped five sets of DCs, including damage to single simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure Specifically, nine vibration MPs were set up at the joint positions, taking into account th shape and characteristics of the structure. From above, the MPs were distributed at th central point of the dome and in the directions of east, west, south, and north. In the ver tical direction, the MPs were evenly distributed at the top, middle, and bottom of th structure. During the time history analysis on the FE models, acceleration response dat were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In th FE analysis process, stiffness reduction was applied to local timber beam and steel cabl elements to simulate damage or stiffness degradation of the components. This paper de veloped five sets of DCs, including damage to single and multiple components. The fiv simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of t Specifically, nine vibration MPs were set up at the joint positions, takin shape and characteristics of the structure. From above, the MPs were central point of the dome and in the directions of east, west, south, an tical direction, the MPs were evenly distributed at the top, middle, structure. During the time history analysis on the FE models, accelera were extracted from the nine MPs. The layout of the MPs was shown FE analysis process, stiffness reduction was applied to local timber be elements to simulate damage or stiffness degradation of the compone veloped five sets of DCs, including damage to single and multiple com simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were Specifically, nine vibration MPs were set up at the jo shape and characteristics of the structure. From abo central point of the dome and in the directions of eas tical direction, the MPs were evenly distributed at structure. During the time history analysis on the FE were extracted from the nine MPs. The layout of the FE analysis process, stiffness reduction was applied elements to simulate damage or stiffness degradatio veloped five sets of DCs, including damage to single simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.  In practical engineering, vibration sensors were used for SHM of the dome structure. Specifically, nine vibration MPs were set up at the joint positions, taking into account the shape and characteristics of the structure. From above, the MPs were distributed at the central point of the dome and in the directions of east, west, south, and north. In the vertical direction, the MPs were evenly distributed at the top, middle, and bottom of the structure. During the time history analysis on the FE models, acceleration response data were extracted from the nine MPs. The layout of the MPs was shown in Figure 6. In the FE analysis process, stiffness reduction was applied to local timber beam and steel cable elements to simulate damage or stiffness degradation of the components. This paper developed five sets of DCs, including damage to single and multiple components. The five simulated DCs are detailed in Table 4.   Table 4. In this paper, the process of establishing a time series model was : damaged by 80%.  Table 4. In this paper, the process of establishing a time series model was demonstrated using the acceleration time history data of MP M, which is located at the top of the dome structure.
Before establishing the ARMA model, it is necessary to standardize the sample dataset. The normal distribution probability plot of the processed data, which was shown in Figure 7c, indicated that the processed sample is a stationary, normal, and zero-mean time series, satisfying the prerequisite for establishing a time series model. After preprocessing the data, the ACF and PACF of the data sample were calculated, as shown in Figure 7a,b. The ACF and PACF gradually decreased to within the 95% confidence interval, and both exhibited tailing characteristics. According to Table 1, an ARMA model could be established for the data sample. Before establishing the ARMA model, it is necessary to standardize the sample dataset. The normal distribution probability plot of the processed data, which was shown in Figure 7c, indicated that the processed sample is a stationary, normal, and zero-mean time series, satisfying the prerequisite for establishing a time series model. After preprocessing the data, the ACF and PACF of the data sample were calculated, as shown in Figure 7a,b. The ACF and PACF gradually decreased to within the 95% confidence interval, and both exhibited tailing characteristics. According to Table 1, an ARMA model could be established for the data sample. Subsequently, the AIC and LAR method were used to determine the order and estimate the parameters of the model. Figure 7d shows the AIC results of the ARMA(2m,2m−1) model for the first 30 orders. Based on the AIC results, an ARMA (20,19) model was established for the MP M, and its estimated parameters are shown in Table 5.  Subsequently, the AIC and LAR method were used to determine the order and estimate the parameters of the model. Figure 7d shows the AIC results of the ARMA(2m,2m−1) model for the first 30 orders. Based on the AIC results, an ARMA (20,19) model was established for the MP M, and its estimated parameters are shown in Table 5.  Table 4. To alleviate the computational burden arising from the high AR order, dimensionality reduction could be achieved by conducting PCA on its parameters. Based on the principle of PCA, the AR parameters of each MP were analyzed for five distinct DCs. Table 6 presents the eigenvalue and ACR obtained through PCA for the data of each MP. The eigenvalue quantifies the extent to which each principal component explains the variation in the data, while the ACR represents the proportion of total variance explained by the first n principal components. An ACR exceeding 85% indicates that the primary components already encompass most of the feature information [26]. Based on the PCA results, the eigenvalues of the first and second order AR coefficients for each MP were nearly 10 and 4, respectively. Furthermore, the ACRs for the first three orders at each MP all exceeded 90%, surpassing the 85% threshold. In some cases, the ACR even reached approximately 99%. Consequently, it can be inferred that employing the first three order AR coefficients to construct DSF can yield sufficiently accurate outcomes in theory. The acceleration time series data from each MP under the healthy state (DC 1) was subjected to windowing processing, resulting in two equally divided sections. The first 60 s of the resulting time series data, comprising 6000 sample data points, was utilized as training sample data, while the last 60 s, also comprising 6000 sample data points, was used as reference sample data. The remaining data obtained from other DCs were unknown-state sample data. The first three order AR coefficients of the ARMA model were extracted for each sample data. A principal component matrix was established and the MD was calculated between the principal component matrices of a single MP under different DCs. Corresponding DSF for damage identification was constructed according to Equation (14). The DSF for each MP is presented in Figure 8. training sample data, while the last 60 s, also comprising 6000 sample data points, was used as reference sample data. The remaining data obtained from other DCs were unknown-state sample data. The first three order AR coefficients of the ARMA model were extracted for each sample data. A principal component matrix was established and the MD was calculated between the principal component matrices of a single MP under different DCs. Corresponding DSF for damage identification was constructed according to Equation (14). The DSF for each MP is presented in Figure 8.
(e) (f)  Based on the results of damage identification, the following observations could be made: (

Damage Identification Based on SHM Data
This paper presented a damage identification study of the dome structure at Taiyuan Botanical Garden, utilizing measured data from vibration sensors. The structure employed glued laminated timber beam as its primary load-bearing component, and exhibited a large span, complex natural vibration characteristics, and susceptibility to wind vibration, earthquakes, and other impacts. To ensure the structure's safety, long-term SHM has been conducted since its service. Based on the structure's characteristics and FE analysis results, nine vibration MPs shown in Figure 6 were selected to monitor the structure's vibration behavior. These MPs were equipped with the 2D001 magnetoelectric vibration sensors produce by Donghua Testing Technology Co. Ltd. (Taizhou, China), with a maximum sampling frequency of 100 Hz. The data acquisition device is shown in Figure 9, and the technical parameters of the 2D001 sensor are shown in Table 7.
Sensors 2023, 23, x FOR PEER REVIEW vibration behavior. These MPs were equipped with the 2D001 magnetoele sensors produce by Donghua Testing Technology Co. Ltd. (Taizhou, China imum sampling frequency of 100 Hz. The data acquisition device is show and the technical parameters of the 2D001 sensor are shown in Table 7.  Under normal working conditions, the only dynamic load on this struct load. However, due to the fact that the dome is covered with glass, the stru posed to sunlight radiation during the day, resulting in temperature gradients that may have an uncertain impact on the sampling data and final analysis fore, this study selected measured data at night as the sample data. Specificall data from the initial stage of structural service (October 2020) were chosen a state sample data, while data after two years of service (October 2022) were unknown-state sample data, with a sample length of 120 s. Hence, ARMA performed and the health status of the structure after two years of service w Figure 10 displays partial data for the initial stage and after two years of serv  Under normal working conditions, the only dynamic load on this structure was wind load. However, due to the fact that the dome is covered with glass, the structure was exposed to sunlight radiation during the day, resulting in temperature gradients on the surface that may have an uncertain impact on the sampling data and final analysis results. Therefore, this study selected measured data at night as the sample data. Specifically, acceleration data from the initial stage of structural service (October 2020) were chosen as the healthy-state sample data, while data after two years of service (October 2022) were selected as the unknown-state sample data, with a sample length of 120 s. Hence, ARMA modeling was performed and the health status of the structure after two years of service was evaluated. Figure 10 displays partial data for the initial stage and after two years of service.
The results of the calculation indicated that the ACF and PACF of the model are tailing, which satisfies the prerequisite for ARMA modeling. The AIC was employed to determine the order of the ARMA(2m,2m−1) model for the measured data. The findings reveal that the ARMA(28,27) model had the lowest AIC value of −1.4235 among the first 30 orders, indicating that this order scheme can achieve a high degree of fitting. The vertical acceleration data of nine vibration sensors at each MP were preprocessed, and ARMA(28,27) models were established. PCA was performed on the AR parameters using a set of measured data from each MP under both healthy and unknown states. The analysis results showed that the first 5 order AR parameters of any ARMA model contained 90% of the model's damage characteristic values. Therefore, the first 5 order AR coefficients were utilized to construct subsequent DSF. The results of the calculation indicated that the ACF and PACF of the model are tailing, which satisfies the prerequisite for ARMA modeling. The AIC was employed to determine the order of the ARMA(2m,2m−1) model for the measured data. The findings reveal that the ARMA(28,27) model had the lowest AIC value of −1.4235 among the first 30 orders, indicating that this order scheme can achieve a high degree of fitting. The vertical acceleration data of nine vibration sensors at each MP were preprocessed, and ARMA(28,27) models were established. PCA was performed on the AR parameters using a set of measured data from each MP under both healthy and unknown states. The analysis results showed that the first 5 order AR parameters of any ARMA model contained 90% of the model's damage characteristic values. Therefore, the first 5 order AR coefficients were utilized to construct subsequent DSF.
The acceleration time series of each MP in a healthy state was extracted and partitioned into two equal segments by windowing. The first 60 seconds' 6000 sample data points were designated as training sample data, while the last 60 seconds' 6000 sample data points were designated as reference sample data. The remaining data were used as unknown-state sample data. The MD was computed to generate DSF for damage identification.
According to the analysis results, the MDs of each sample were not significantly different, and the DSF of all MPs were close to 1.0, indicating that the structure was almost not damaged after two years of service. The results of MD and DSF of each MP were shown in Figure 11. In addition, it could be found that the DSF of MP S2 was 1.097, while the maximum value of DSF of other MPs was 1.062, and the average value was 1.046, indicating that the damage degree of the component at S2 was slightly higher than that of other MPs. Figure 12 shows the static FE analysis results of the dome. Under the combined action of dead load and roof live load, the stress level of the timber components around S2 was higher than that in other areas. Moreover, real-time dynamic loads such as wind loads and temperature effects might have adverse effects on local areas. These adverse effects might be the potential factors leading to the increase of DSF at S2. Therefore, it is necessary to pay attention to the SHM data of this dome structure, and focus on the acceleration data of MP S2, so as to timely warn when the structure suffers serious damage. The acceleration time series of each MP in a healthy state was extracted and partitioned into two equal segments by windowing. The first 60 seconds' 6000 sample data points were designated as training sample data, while the last 60 seconds' 6000 sample data points were designated as reference sample data. The remaining data were used as unknown-state sample data. The MD was computed to generate DSF for damage identification.
According to the analysis results, the MDs of each sample were not significantly different, and the DSF of all MPs were close to 1.0, indicating that the structure was almost not damaged after two years of service. The results of MD and DSF of each MP were shown in Figure 11. In addition, it could be found that the DSF of MP S2 was 1.097, while the maximum value of DSF of other MPs was 1.062, and the average value was 1.046, indicating that the damage degree of the component at S2 was slightly higher than that of other MPs. Figure 12 shows the static FE analysis results of the dome. Under the combined action of dead load and roof live load, the stress level of the timber components around S2 was higher than that in other areas. Moreover, real-time dynamic loads such as wind loads and temperature effects might have adverse effects on local areas. These adverse effects might be the potential factors leading to the increase of DSF at S2. Therefore, it is necessary to pay attention to the SHM data of this dome structure, and focus on the acceleration data of MP S2, so as to timely warn when the structure suffers serious damage.

Conclusions
Structural damage identification is a crucial analysis method in SHM, with significan scientific and engineering applications. To address the challenges of incomplete moda information, high modal density, and large degrees of freedom in large-span spatial lattice structures, this paper proposed a structural damage identification method based on the ARMA model and PCA method for structures of this kind. Based on the above research the following conclusions can be drawn: (1) This article was based on a time series model and introduced PCA for extracting principal components. In this study, only the first three AR coefficients were required to achieve damage identification and it demonstrated high accuracy. (2) The MD is advantageous compared to other distance methods in constructing DSF because it considers the differences and correlations in the variability of variables fo each observation. Additionally, it can effectively calculate the influence of scale on different measurement values. (3) Verified by the FE model, the damage identification method proposed in this article was found to be sensitive to structural damage. It could, to a certain extent, accurately

Conclusions
Structural damage identification is a crucial analysis method in SHM, with significant scientific and engineering applications. To address the challenges of incomplete modal information, high modal density, and large degrees of freedom in large-span spatial lattice structures, this paper proposed a structural damage identification method based on the ARMA model and PCA method for structures of this kind. Based on the above research, the following conclusions can be drawn: (1) This article was based on a time series model and introduced PCA for extracting principal components. In this study, only the first three AR coefficients were required to achieve damage identification and it demonstrated high accuracy. (2) The MD is advantageous compared to other distance methods in constructing DSF because it considers the differences and correlations in the variability of variables for each observation. Additionally, it can effectively calculate the influence of scale on different measurement values. (3) Verified by the FE model, the damage identification method proposed in this article was found to be sensitive to structural damage. It could, to a certain extent, accurately locate the site of structural damage and reflect the extent of the damage. Therefore, it can be effectively utilized for practical structural damage identification. (4) Based on SHM data, the paper identified the damage of the structure after two years of service. The results indicated that the structure is in a relatively healthy state, but the DSF of MP S2 is slightly higher than those of other MPs. In the future, it is necessary to pay attention and give timely warnings if necessary.
However, there are still some shortcomings in this article. The next step should be to conduct more in-depth research on the rationality of the number and distribution of measurement points, as well as the refinement of the FE model correction.

Data Availability Statement:
The data presented in this study are available from the first author, (wangguoqing0066@link.tyut.edu.cn), upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.